Modeling Longitudinal Oscillations of Bunched Beams in Synchrotrons 
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Longitudinal oscillations of bunched beams in synchrotrons have been analyzed by accelerator 
physicists for decades, and a closed theory is well-known The first modes of oscillation are the 
coherent dipole mode, quadrupole mode, and sextupole mode. Of course, these modes of oscillation 
are included in the general theory, but for developing RF control systems, it is useful to work with 
simplified models. Therefore, several specific models are analyzed in the paper at hand. They 
are useful for the design of closed-loop control systems in order to reach an optimum performance 
with respect to damping the different modes of oscillation. This is shown by the comparison of 
measurement and simulation results for a specific closed-loop control system. 

PACS numbers: 29.20.dk, 29.27.-a 



I. INTRODUCTION 

According to the standard theory, longitudinal bunch 
oscillations are characterized by different mode numbers. 
The mode number m £ {1,2,3,...} describes the shape 
of an individual bunch in phase space whereas the mode 
number n € {0, 1, M — 1} describes the phase relation 
between the oscillation of M individual bunches in case 
there is coupling between the bunches (e.g. n — speci- 
fies the case that all bunches are oscillating in-phase). 

In the paper at hand, only single-bunch oscillations are 
investigated. 

In longitudinal phase space, the within-bunch modes 
m may be visualized by a symmetrical polygone with 
m rounded corners. This symmetry of course requires 
appropriate scaling of the phase space coordinates. 

If the bunch is located in the linear region of the 
bucket, i.e. if the revolution frequency in phase space ap- 
proximately equals the synchrotron frequency fs = 1/ Ts 
in the bucket center, it is obvious that the phase space 
distribution will be the same one as the initial one after 
the time Ts/m. Therefore, it is clear that the beam signal 
(projection of phase space ensemble onto the time axis) 
will contain spectral lines at mfs- Since the same bunch 
re-appears in a ring accelerator after the revolution time 
Tr = 1//r, spectral lines are observed at Q 

f= P Mf R ±mf s for p£ {0,1,2,...}. 

Therefore, the existence of a spectral line with a specific 
value m may be used to define the mode of oscillation. 
Using this definition, one concludes that a quadrupole 
mode is present if a spectral line fitting to m = 2 is 
observed. 

Considering only single-bunch oscillations is not the 
only simplification that is made in this paper. Many 
other effects are not taken into account which are essen- 
tial for high-current beam acceleration: 



• No real beam instabilities are considered. This 
means that no beam impedances are taken into ac- 
count, and no beam loading is present. 

• No space charge effects are considered. 

• No coupling between longitudinal and transverse 
beam dynamics is assumed. 

The exclusion of these effects is usually not relevant for 
the first design of a closed-loop control system. Another 
reason for these simplifications is that in spite of them 
some effects occur which deviate from explanations that 
can be found in literature. Therefore, it is clear that one 
has to be even more careful when observations are gener- 
alized by adding further physical phenomena of practical 
importance. 

In the following sections, several models for longitu- 
dinal single-bunch oscillations will be presented that de- 
scribe the phenomenon from different points of view. Af- 
terwards, a specific closed-loop control system for damp- 
ing both, dipole and quadrupole oscillations is presented. 
For this system, measurement results are compared with 
simulation results. 



II. MODELS 

In the following, different models for describing longi- 
tudinal bunch oscillations are presented. 



A. Mode Definition by Bunch Shape in Phase 
Space 

Assuming that the phase space is scaled in such a way 
that a matched bunch in the linear region of the bucket is 
a circle, the m-th mode may be described by the formula 



r = 1 + €b sin(m ip). 



(1) 
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Here, (r, tp) are the polar coordinates of the bunch con- 
tour. It is obvious that an unperturbed bunch is obtained 
for m = 0. 

Fig. [T] shows the bunch contours for €b = 0.1. 
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FIG. 1. Phase space contour of longitudinal oscillation modes 
B. Simplified Linear Differential Equations 



We assume that the RF voltage is modulated according 



to 



where 



i(t) = fio(t)(l + e(t)) sin(ip(t) - Ap gap (t)) (2) 



URF(t) dt. 



By definition, the reference particle arrives at the ac- 
celerating gap when ip(t) = <pii(t) + 2nk is valid (the 
integer k denotes the bunch repetition number); for a 
non-synchronous particle, the arrival time is defined by 
ip(t) = <pR\(t) + 2nk + Aip(t). The magnetic field B in the 
bending dipoles and the quantities uq, uirf and ipR are 
chosen in such a way that the reference particle follows 
the reference path. All these quantities vary slowly with 
time in comparison with the synchrotron oscillation. The 
modulation functions e(t) and Ap gap (t), however, may 
vary faster. 

Thus, the nonlinear differential equations are 



AW 



2TThf] R 

T R p%w R 



AW 



Quo 



Tr 

sin ipR] 



[(1 + e) • sin(ip R + Aip - Aip gap ) 



(3) 



(4) 



Here, h denotes the harmonic number, Tr the revolution 
time, Wr — 7tioCq7_r the total energy, and /3r, jr are 
the relativistic Lorentz factors of the reference particle. 
With the transition gamma jt, VR ~ V7t — V7ii ^ s * ne 
phase slip factor. Q is the charge of one single particle. 
AW and Aip are the energy and phase deviations of a 
non-synchronous particle with respect to the synchronous 
reference particle (cf. 0). 

For small values of \Aip — Aip gap \ < 1 we have 

aw ~ Quocos <p R ( etanipR 
Tr ^ £ > I ^ ~ lfaap 1 + e 



Defining 



A<Pgap = A<^ 



gap 



1 + 6 



tan ipR 



leads to 



AW 



Juqcos pr 
Tr 



(l + e)(Atp- Aipg ap ). 



(5) 



(6) 



By a combination of equations ([3]) and (O, we get 
2i:hriRQuQCOs pr. 



Aip = 



TrPrWr 

The synchrotron frequency is defined by 



(l + e)(Ap- Ap gap ). 



which yields 



< 2whQu (-r)R cos <p R ) 



WrWr 



(7) 



Aip + u s (l + e)Ap = ui s (l + e)Aip gap . 
By using the new variables 

x — Aip y = CAip 

we find: 

± = h v 

y = + e)(x - Ap gap ) 

For vanishing excitations with e = and Aip gap = 0, we 
now require the trajectories to be circles: 

x = cos(ujst) => y — Cx — —Cu>ssin(ujst) =>• C = 

LUs 

Thus, we obtain (note that C and vary slowly - 
therefore we neglect the time derivative): 

x = -usy (8) 

y = cJs(l + t)(x ~ Aipgap) (9) 

C. Behavior of Particle Bunches 

Whereas equations ([8|) and ((9]) are valid for individual 
particles, we now consider bunches with N particles. 

1. Phase Oscillations 



For the mean values, we find: 



1 N 

-Y 



Xk y 



fc=i 



l N 



k=l 



1 N 1 N 



fc=i 



fc=i 
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1 N 1 N 

^^2ilk^ u s {l + e)—^2{x k - Aip gap ) 



k=l 



fc=l 



y = us(l + e)(x - Aip gap ) (11) 



x = -Lj s y = + e)(ar - Aip gap ) (12) 

This equation for the bunch center has the same form as 
equation {3 for the individual particles. 



2. Amplitude Oscillations 
We define the following quantities: 

N , JV 



1 N 1 W 1 * 



fc=i 



fc=i 



fc=i 



j JV 1 AT 

«x = -^(ii-if = — ^(x 2 k -2x k x + x 2 ) = a x -x 



k=l 



k=l 



v y = a y - y 

Please note that v x corresponds to the variance of the 
quantities x k if a division by N — 1 is used instead of 
the division by N . Since we are only interested in large 
numbers N, this difference is negligible. 

The quantity ^JvZ represents the bunch length whereas 
represents the height of the bunch (this will later be 
analyzed in detail). We get: 

1 N 

a x = — 2J 2xfeifc = ~2uj s £, 

k=l 



1 N 1 N 

>',: = —>2ykVk = 2o;g(l + e)—y]yk(xk - Afi gap ) 



k=l 



k=l 



d y = 2ws(l + e)(£ - yAif ga p) 



1 N 

i = j^^2(±kyk + x k y k ) 



k=l 



= -bj s a y + uj s {1 + e)(a x ~ xA(p gap ) 

Vx = d, x — 2xx = —2ugS, + 'Zusxy — —2ioga (13) 

Here we defined a — £ — xy in order to have the same 
form in the expressions for a x and for v x . We find: 

u y = ay - 2yy = 2u s (l + e)(£ - yAtp gap ) - 
- 2w s (l + e)y(x - A(p gap ) 



=> v y = 2ws(l + e)a 

a = £ — xy — xy = 

= -Uga y + w s (l + e)(a x ~ xACp gap ) + 
+ oj s y 2 - w s (l + e)x(x - Aip gap ) 

=> a = -usVy +uj s {1 + e)v x 



(14) 



(15) 



Now we are able to derive a differential equation for v x , 
i.e. for the bunch length oscillation. 

Combining equations (fi"3|) and (fT4l) yields: 

v v = -(l + e)v x (16) 

We now combine equation (fT3f with equation f| 15|) : 

^ = 2ugV y - 2w|(l + e)uj (17) 

=>'v x = 2u!gi>y — 2cjg(l + e)^ — 2u> s cv x 
Using equation ([T^]). we finally get: 

v x = -4cj|(1 + e)uai - 2u>gev x 



(18) 



Please note that for e = 0, the standard differential equa- 
tion 



v x 4- (2ujg) 2 v x — const. 



(19) 



is obtained which corresponds to an oscillation with 
the frequency 2tog. Due to the linearization, an initial 
quadrupole oscillation will continue forever. 

Now we derive the differential equation for v v , i.e. for 
the amplitude oscillation. Equation (fT4)) yields: 

= 2u)gu 



1 

The time derivative is 

%(! + e) ~ 
(l + e)2 



-2u s v y + 2cj§(1 + e)v x 



where we used eqn. (jT5j) on the right side. We divide by 



(1 + 6) 2 (l + e)3 



-2ui 



1 



2lo 2 s v x 



Now another time derivative leads to v x on the right side 
such that we can use eqn. (|16[) to eliminate v x completely. 
After some steps, one obtains: 



3wj,e 



2ujgCV y 



V y \ AljJg{l + < ) 



(20) 



This differential equation for v y differs from eqn. (|18[) 
for v x only by terms that are of higher-order with re- 
spect to e. Furthermore, the sign of the excitation term 
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2uj 2 s ivy is different for v x and v y which matches the ex- 
pectations since the bunch is short when its amplitude 
is high whereas the bunch is long when its amplitude is 
small. 

Please note that we have not introduced any approxi- 
mations to derive the differential equations (fl8l) and (|20|) 
from equations © and ([§]). 

According to [4(, these differential equations have the 
following solution: 

v x = C x iw 2 xl + C x2 w xl w x2 + C x3 wl 2 (21) 
v y = C y iw 2 yl + C y2 w yl w y 2 + C y3 w 2 2 (22) 

The functions w x i and w x2 are the linearly independent 
solutions of 

w x + + e)w x = 0, 

whereas the functions w y \ and w y2 are the linearly inde- 
pendent solutions of 

e 2 
w y - Y^r^y + u s( l + € ) w y = °- 

In the trivial case e = 0, we may choose 

w x i = w y i = cos(uJst), w x2 = w y2 = sin(bj s t) 

as a solution. Due to equations (f2~Tj) and ([2"2")) , v x and v y 
will oscillate with twice the frequency which is in com- 
pliance with eqn. ([T9| . 

In appendix [K\ it is shown that the following differ- 
ential equations are valid as an approximation for small 
deviations from the matched bunch shape: 

x + uigx = uigAtpgap (23) 
v x + Aujg(v x - v ) = -2u> 2 s v e. (24) 

D. Revolution Time in Phase Space 

The revolution frequency of off-center particles in the 
nonlinear bucket in phase space is given by 

7T 

fs .nonlinear (A<p) = fs t (25) 

2x(sin^) 

(cf. 3). Here, A<p is the maximum phase deviation of 
the particle and K is the complete elliptic integral of the 
first kind. Due to the longitudinal emittance of the bunch 
in the nonlinear bucket, the quantities x and v x will not 
oscillate with the frequency fs and 2/5, respectively, but 
with reduced frequencies fs,eff and 2fs, e ff- 

We now analyzed how fs,eff depends on the size of 
the bunch. For this purpose, the solutions of the lin- 
earized ODEs (j2"3"]l and (|24p were compared with a non- 
linear tracking simulation for an initial bunch with an 
example shape given by 

r = r [1 + 0.1 cos(^) + 0.2 cos(2^)] . 
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FIG. 2. Effective phase deviation A<p e ff for different bunch 
sizes and different modes of oscillation ( X: dipole mode, 
o: quadrupole mode, black: homogeneous distribution, grey: 
Gaussian distribution) 

The longitudinal emittance is determined by the parame- 
ter ro . Due to the mismatch of the bunch, v x and v y per- 
form oscillations around the average values v x q and v y o, 
respectively. Please note that one has v x q < v y o since 
the nonlinear bucket has the shape of an eye instead of a 
circle. For each run of the simulation, the average value 
v x q was determined as a measure for the bunch length. 
Each run of the nonlinear tracking simulation also leads 
to a certain oscillation frequency, and we determined Aip 
such that this frequency matches fs, nonlinear (Aip) for x 
or 2 f s , nonlinear (Aip) for v x , respectively. As a result, the 
effective phase deviation Aip =: Ap e ff and v x q are de- 
termined for each simulation run. 

Fig. [2] shows the relation between these two quanti- 
ties. It is obvious that the effective synchrotron frequency 
fs,eff is approximately determined by 

fs,eff — fs, nonlinear {Aipeff) (26) 

if 

Aip e ff = 2y/v^ (27) 

is used. 

E. Interpretation of x and v x 

According to the equations 

x w —<pi, v x w 2 In— — 2 ln—^- (28) 

the mean value x and the variance v x may approximately 
be converted into the phase ip\ and the amplitude A\ of 
the fundamental harmonic component. Aq is the zeroth 
Fourier component, and i is the DC component of the 
beam signal. The approximation (|28p was derived for an 
elliptical Gaussian bunch with a large number of parti- 
cles. It also holds for a homogeneous distribution. 
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TABLE I. Simulation parameters 



216.72 m 
5.45 

4„ Ar 18+ 

11.4 MeV/u 
10 kV 



Synchrotron circumference Ir 
Transition gamma 7t 
Ion species 
Kinetic energy 
RF amplitude uo 
Harmonic number h 8 
Synchrotron frequency fs 3312 Hz 

Revolution time Tr 4.66 /is 

Parameters resulting from the initial conditions 
v x o 0.773 
A<p eff = 2v^ 1.7584 
fs,eff/fs 0.8087 




100 
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50 
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100 



FIG. 3. Mean value x and variance v x versus time for 
ws,eff = 0.8 uus- The initial bunch is a homogeneous dis- 
tribution according to r = 1.6 [1 + 0.1 cos(^j) + 0.2 cos(2</p)] . 
No excitation is applied ( x : nonlinear tracking simulation, 
+: result of the FFT using eqn. {35]), A: nonlinear ODE, V: 
linearized ODE). 



We will show now that equation (|28[) is also a good 
approximation for larger bunches and even if a significant 
amount of filamentation is present. 

For the parameters shown in Table U solutions for the 
following models were generated numerically: 

1. A nonlinear particle tracking simulation. The pa- 
rameters x and v x are calculated. 

2. An FFT analysis of the previous result was made 
which leads to the phase ipi, the amplitude A\ and 
the DC component Ao/2 = i. Based on these pa- 
rameters, an approximation for x and v x is calcu- 
lated using eqn. ([28]). 

3. Solution of the nonlinear ordinary differential equa- 
tions (ODE) (H2) and (HU). 

4. Solution of the linearized ODE (E! and (l24l). 



For the last two models 3 and 4, the effective syn- 
chrotron frequency fs.eff was used instead of /<?. 

Fig. H shows the result for a mismatched bunch without 
additional excitation, and Fig.[5]shows the corresponding 
phase space plots. Fig. 2] shows the result for the same 
mismatched bunch with an additional excitation. These 
examples show the following: 

• The formula (|2"8"|) describes the oscillation very ac- 
curately, but there is a small DC offset which may 
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FIG. 4. Mean value x and variance v x versus time for 
ws,eff = 0.8 ljs- The initial bunch is a homogeneous distri- 
bution according to r = 1.6 [1 + 0.1 cos(^) + 0.2 cos(2y>)]. An 
excitation e = 0.04 sin(2ojs. e //t + 7r) is applied ( X: nonlinear 
tracking simulation, +: result of the FFT using eqn. (|28[l . A: 
nonlinear ODE, V: linearized ODE). 
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FIG. 5. Phase space distribution at the beginning (left di- 
agram) and at the end (right diagram) of the simulation. 
The initial bunch is a homogeneous distribution according 
to r — 1.6 [1 + 0.1 cos(y) + 0.2 cos(2</?)] (mismatch of both, 
bunch center and bunch length). No excitation is applied 
(same case as in Fig. [3j). 



be relevant for small oscillation amplitudes. As we 
will see later, the DC component is usually not of 
interest for the control loop design. 

• The models 1 (marker x in the diagrams) and 
2 (marker + in the diagrams) match very well 
even for comparatively large bunches and for mis- 
matches including filamentation. This clearly 
shows that the quantities x and v x may be used 
instead of (p% and A\ if the approximation (|28[) is 
used. 

• The models 1 and 2 include Landau damping since 
they are based on nonlinear tracking equations. 
The models 3 (marker A) and 4 (marker V) can- 
not show Landau damping since they are based on 
linear tracking equations. 

• During the first oscillation period, all models match 
very well which indicates that all models may be 
used for designing feedback systems. 

• Fig. U shows that the excitation with 2fs, e ff ini- 
tially leads to a damping of the amplitude oscilla- 
tion. The initial damping rates of all four models 
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are similar. This is a further confirmation that the 
models may be used for control loop design. 



F. Phase and Amplitude Oscillations of the 
Normal Modes 

We now consider the normal modes introduced in sec- 
tion III Al If one assumes a homogenous particle distri- 
bution inside the bunch contour defined by equation ([1} , 
one may calculate the phase and amplitude oscillations 
analytically. The integrals that have to be solved are the 
following ones: 

Mo = I / r dr dip 
Jo Jo 



x = / / x r dr dip 

Mo Jo Jo 



2tt rTmax 

'i — yt~ I / y r dr d v 

Mo Jo Jo 




3 4 5 6 7 8 

Time/\is 

FIG. 6. Visualization of the partial sums of the series in 
eqn. ^ for Tr = 1 fis, T s = 8 (is, At = 0.3 fis (solid 
line: summation up to n = 20, dashed line: summation up to 
n = 50) 
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y 2 r dr dip, v y = a y - y 2 



o Jo 

Here, the following equations were used: 

r ma x = l + e B sin [m (ip - ip )} 

x — r cos ip y — r sin ip 
The following results are obtained: 

for m = 



Mn = 



7T 
7T fl 



if 2 

2 B 



for m > 



FIG. 7. Dipole oscillation of cos 2 -shaped bunches for Tr = 
1 /xs, Ts = 8 fxs, At — 0.3 fis, f2 = 18/Tr (summation up to 
n = 100) 



modulation with the frequency 2fs (since v y is periodic 
with respect to 2ipo), but this is a parasitic effect caused 
by the fact that the bunch shape is only approximately 
circular. The slight deformation causes the modulation 
of the order e 2 B . Apart from this exception, only the 
mode m — 2 shows an amplitude modulation with the 
frequency 2f$ of the order es- The modes m > 2 neither 
show a phase modulation nor an amplitude modulation. 



Vy = < 



1 t o e ? 



for m = 1 
for m ^ 1 



for m = 1 



(eb+|<:b) cos vo 
for m^l 



for 







3,2 i 1A _1_ 1 ,6 



1 tTj'jT, 
4 



*B + T 



e„-o£ n cos 



(2yo) 



for m = 1 



4 

_L 

L 4 



(1+14) 

1 l+34 + §4 + (2 eB + §4)sm(2y ) for ^ = 2 



2 fc B 



for m > 2 



We see that only the mode m — 1 shows a phase mod- 
ulation with the frequency fs (since x is periodic with 
respect to ipo). The mode m = 1 also shows an amplitude 



G. Spectrum of Dipole Oscillation 

In appendix [B] we prove the formula 



Sit — kTft — At sin 2nk— I 

k— — OG ^ \ / 



— + ^2 A n cos{nuj s t) (29) 



Tr 



for At < T R /2 where 



1 + cos [ riajg 

k=l 



fcTfj + At sin 



. 2nkT R 



Ts 



and Ts/Tr is a positive integer. 

Fig. [B] shows the graph of two different partial sums 
of the series on the right side of eqn. (|2"9"|) for one period 
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Ts = 2tt/ujs- It is obvious that the peaks become the 
higher the more terms are added. One also sees that 
the pulse density is high in the middle whereas it is low 
at the beginning and at the end of the period. This 
is the expected behavior for the dipole oscillation under 
consideration. 

A real beam signal consists of pulses with finite height 
and length instead of the Dirac pulses. Let us assume 
that such a single finite pulse centered at t = is given 
by the function x(t). If we write the above-mentioned 
train of Dirac pulses as 

+00 

h(t)= y, Kt-n), 

we get the convolution 

y(t) = h(t) * x(t) = X (t- T k), (30) 



fe=-c 



+00 



k=— 00 



which corresponds to the desired beam signal. This signal 
still shows dipole oscillations if Tfc is implicitly defined by 
eqn. Since 



Ht)= J2 



c (h) e jnu:st 



n— — 00 



is given by a Fourier series, the corresponding Fourier 
transform is 



+00 



n— — 00 



where = 4f-. The Fourier transform of eqn. (1301) is 

Y(juj) = H(ju)X(ju) = 

+ 00 

= 2ir Y c^X{jnu, s )S(cj-nu; s ). 



n— — 00 



This is again a Fourier series whose Fourier coefficients 
are given by 

For beam pulses that have the form 



x(t) = \ C ° S v , 



2 1 Ut) for - 7T < fit < 7T 



elsewhere 



one may — as a simple but still realistic example 
derive the Fourier transform 



HT^ft for M + n 



2Q 



for \cj\ = f2 



An example for such a bunch train with dipole oscilla- 
tions is shown in Fig. [7J In this example, we have chosen 



an unrealistic synchrotron period of Ts — 8 fis in order 
to be able to show the result in a diagram. For a re- 
alistic synchrotron period of Ts = 1000 fis, the Fourier 
coefficients are shown in Table [TTJ Please note that the 
following condition has to be fulfilled since the bunch has 
to fit into the bucket even if the maximum bunch offset 
At occurs: 



The following points are observed: 



• The dominant spectral lines do not occur at pfa ± 
fs in general. In our example, this is only valid for 
p=l. 

• Although an ideal dipole mode was modeled, spec- 
tral lines also exist at pfa ± 2f$, p/r ± 3/s, etc. 
This is normally interpreted as a quadrupole, sex- 
tupole, etc. component in the oscillation. 

• The sidebands are not symmetric with respect to 
the harmonics of the revolution frequency. 

When At is reduced significantly in this example (i.e. 
for smaller dipole oscillations), the dominant spectral 
lines are located at pfn ± fs as expected, and the sym- 
metry is also improved. The pure existence of Fourier 
components at pfu ± 2fs, pJr ± 3/s, etc. remains, how- 
ever (their magnitude is smaller than that at pfn ± fs). 

It has to be emphasized that the validity of these ob- 
servations is very general. Even though the example was 
a special one, the formulas for the Fourier series are not 
based on any approximations. 



III. SEXTUPOLE MODE GENERATION 

In this section, we try to excite a sextupole oscillation 
by a phase modulation with 



A(fg 



18 



sin(3ujst) 



which corresponds to a phase deviation of ±10°. The 
parameters defined in Table U are used, and a matched 
bunch (r = 1) is assumed at the beginning. 

After a time t = 1.12 Ts, the phase space distribution 
shown in Fig. [8] is obtained. The left diagram shows that 
no typical sextupole distribution is obtained if linearized 
tracking equations are used. In the right diagram one can 
see that the original nonlinear tracking equations lead to 
a sextupole distribution whose contour is compliant with 
equation ([T]) for m = 3. Matching the contour to the 
particle cloud leads to 

r = 1 + 0.1 sin(tp — n) + 0.1 sin(3ip). 

This shows that both, a dipole and a sextupole mode is 
excited. It has to be emphasized that an excitation of the 
sextupole mode is only possible in the nonlinear bucket. 
This was also verified using the program package ESME 
Q and can furthermore be shown analytically 0. 
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TABLE II. Fourier coefficients oi = 2c„ of a bunch train with dipoie ocsciilations {Tr = 1 /is, Ts = 1000 /is, At = 0.3 /is, 
n = 18/Tfl) 





/r 


2/r 


3/fl 


4/h 


-5/s 


0.0016819 


0.0267036 


0.0545552 


0.0222028 




0.0087296 


0.0618662 


0.0634047 


0.0012028 


-3/s 


0.0355937 


0.1048087 


0.0351619 


-0.0208789 


-2/s 


0.1050056 


0.1050122 


-0.0260749 


-0.0177420 


-1/5 


0.1875566 


0.0065220 


-0.0534534 


0.0114866 


+0/s 


0.0937025 


-0.1014007 


0.0073288 


0.0207027 


+l/s 


-0.187475 


-0.0061292 


0.0533306 


-0.0115380 


+2/s 


0.1055102 


0.1045001 


-0.0263491 


-0.0175559 


+3/s 


-0.0361045 


-0.1048387 


-0.0345397 


0.0208711 


+4/s 


0.0089756 


0.0622868 


0.0629741 


0.0008722 


+5/s 


-0.0017600 


-0.0271106 


-0.0545328 


-0.0217486 



1 

0.5 

» 
-0.5 

-1 




-0.5 0.5 1 

X 




FIG. 8. Phase space at t = l.YlTs for an excitation with 
3/s (Left diagram: linear tracking equations, right diagram: 
nonlinear tracking equations) 



IV. APPLICATION OF MODELS 



Central Control System 



Target 
Frequency 



Target 
Frequency 



I 

6 



Target 
Amplitude 



DSP System A 
Filter 



Beam 
Signal 



Integrator 




Amplitude 


& Filter 


4 


Detection 



Using the models described in section [TH an example 
for a specific longitudinal damping system is presented in 
the following. The theoretical results are compared with 
measurement data. 

A. Beam Phase Control and Quadrupole Damping 
System 

The models presented before were used to analyze an 
RF system that is displayed in Fig. [9j A DDS module 
provides an RF master signal with the frequency pro- 
vided by the central control system. An amplitude con- 
trol loop which is not considered in the paper at hand 
makes sure that the detected amplitude of the cavity RF 
signal matches the target amplitude provided by the cen- 
tral control system. 

A phase detection algorithm in DSP system A com- 
pares the phase of the cavity voltage with the phase of 
the beam signal. This phase difference is then processed 
by a digital bandpass filter with fs belonging to the pass- 
band. After applying a proportional gain, the filter out- 
put is used to modify the DDS frequency. For analyzing 
the loop characteristics, the DDS may be regarded as an 
integrator with respect to the phase. The loop A en- 
sures that coherent dipoie oscillations with a frequency 



FIG. 9. Block diagram of the beam phase control system and 
the quadrupole damping system 



approximately equal to fs will be damped. 

In DSP system B, an amplitude detection algorithm is 
applied to the beam signal. This signal is fed to a digital 
filter with a passband frequency of about 2/s- In contrast 
to the dipoie oscillation damping loop A mentioned be- 
fore, there is no intrinsic integrator in loop B. Therefore, 
an additional integrator is implemented in DSP system 
B. The resulting signal is used to modify the target am- 
plitude in order to damp quadrupole oscillations. 

Both DSP systems need a reference signal which allows 
them to detect RF signals at the relevant frequency. For 
this purpose a reference DDS is used. 



B. Analytic Model for the Damping System 

The first natural step to analyze a control loop sys- 
tem like that shown in Fig. [5] is to linearize the building 
blocks in order to get transfer functions in the Laplace 
domain. For example, eqn. ([2U]) directly leads to the 
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TABLE III. Beam experiment parameters 



1 2 3 

fpass/^fs 



FIG. 10. Stability diagram for damping dipole (left diagram) 
and quadrupole modes (right diagram), Td = 10 fjs, fs = 
3312 Hz 



well-known [7j beam transfer function 



0,(8) = 



Aip B 



gap 



where Aips = x is the beam phase. Eqn. (|24p is more dif- 
ficult to interpret. In a well-designed quadrupole damp- 
ing system, no significant emittance blow-up will occur. 
Therefore, the variance vo is a constant which is defined 
by the longitudinal emittance of the bunches. Hence, by 
defining 

Av x = v x - v 
we may rewrite eqn. (|24[) in the form 

Av x + 4u>g Av x = —2ujgVo£- 
which corresponds to the transfer function 



G 2 (s) 



Av x 
-2v e 



■4 W f 



which is also known from literature [7|. 

By using the transfer function Gi(s), the beam-phase 
control loop was analyzed in 8]. In an analog way we 
took G2(s) to analyze the quadrupole damping loop. In 
both cases, the same type of filter specified in @ was 
used. Fig. [10] shows the region of stability for both, the 
dipole and the quadrupole damping system. Formulas 
for the first one can be found in Q, the stability region 
for the latter one is given by 



Kr 



min K. 
p 



p,crit 



where 



K, 



p,crit j^g £]1 £ 



cosf 



l+ 4 r"/ s for C < 1 and p e {1,3,5,. . .} 



Synchrotron circumference Ir 

Transition gamma 7t 

Ion species 

Kinetic energy 

DC beam current i 

RF amplitude uq 

Harmonic number h 

Synchrotron frequency fs at 10 kV 

Revolution time Tr 



216.72 m 
5.45 

40 Ar 18+ 

11.4 MeV/u 
2 mA 

step from 5 to 10 kV 
8 

3312 Hz 
4.66 us 



TABLE IV. Control loop parameters for measurement and 
simulation 

Delay time dipole damping system 10 /xs 

Delay time quadrupole damping system 10 /is 

Filter frequency quadrupole damping system 9000 Hz 

Filter frequency dipole damping system 3500 Hz 

Sampling frequency of FIR filters 375.44 kHz 



l+trlfs for 1 and p e {1,3,5,. . .} 



and 



J pass 

These results were derived in the same way as described 
in The stability diagram is in compliance with the 
measurement results. 



C. Measurement Results 

In order to verify the models presented here, a closed- 
loop analysis was performed, and the results were com- 
pared with the measurement results obtained in a beam 
experiment. Table Hill shows the beam experiment condi- 
tions. In this experiment, a matched beam in a stationary 
bucket generated by an RF voltage of 5 kV was present. 
By doubling the voltage instantaneously, quadrupole os- 
cillations were excited intentionally. Excitations of this 
magnitude will never occur in practice — they were only 
used to show the functionality of the system and the va- 
lidity of the theory. 

The upper diagram in Fig. QT] shows the measured am- 
plitude of the beam signal. The lowest trace is obtained if 
neither the beam phase control nor the quadrupole damp- 
ing system are switched on. The oscillation frequency is 
lower than 2 fs (period of 200 fis instead of 150 /xs) 
since the nonlinearity of the bucket cannot be neglected 
for the momentum spread of the bunch. Due to Landau 
damping, the oscillation becomes weaker with time. The 
trace in the middle of Fig. Qj] shows that the damping 
time is reduced significantly in comparison with Landau 
damping if the quadrupole damping system is switched 
on. For the last uppermost trace in Fig. Illl not only this 
quadrupole damping system was switched on, but also 
the beam phase control system which damps the coher- 
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jjg -0.2 -0.1 



FIG. 11. Measurement results (upper diagram) and simula- 
tion results (lower diagram) for beam signal amplitude (Up- 
permost trace: both control loops for quadrupole damping 
and for dipole damping on, trace in the middle: only control 
loop for quadrupole damping on, lowest trace: both control 
loops off) 



2000 




FIG. 12. Measurement results (upper diagram) and simula- 
tion results (lower diagram) for beam signal phase (Upper- 
most trace: both control loops for quadrupole damping and 
for dipole damping on, trace in the middle: only control loop 
for quadrupole damping on, lowest trace: both control loops 
off) 



ent dipole mode. It is obvious that both control loops 
work together without negative influence. 

The upper diagram in Fig. fTS] shows the measured 
phase of the beam signal. A phase offset was intro- 
duced for each trace in order to have the same order as in 
Fig. [TTJ Therefore, the lowest trace again corresponds to 
both control loops being off which leads to Landau damp- 
ing only. As expected, the oscillation period of 400 /is 
equals twice the period of the quadrupole oscillation. The 
uppermost trace shows that the coherent dipole oscilla- 
tion is damped faster than the Landau damping time if 
both control loops are switched on. The trace in the mid- 
dle indicates that switching on the quadrupole damping 
system while disabling the dipole damping system in- 
creases the excitation of the coherent dipole mode (in 
comparison with the case that both loops are disabled). 




t/,ttS 



FIG. 13. Measured beam signal (Dotted line: before voltage 
jump, solid line: first maximum after voltage jump) 



D. Simulation 

The experiments were compared with nonlinear track- 
ing simulations. The simulation program used nonlinear 
discrete mapping equations in Atp and AW/ujrf for the 



longitudinal dynamics. The particle positions in phase 
space were then converted to the (x, y) plane with the 
variables x — Aip and y = —Atp/us- The beam cur- 
rent signal was calculated as a histogram using bins on 
the A(/3-axis. The beam signal amplitude and phase were 
obtained by an FFT of the beam current signal. 
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TABLE V. Variances of measured beam before voltage step 





V x = V x 


1),. — 7J»,n 


Vn — „ / V ,-V>, 

^ U \/ « J/ w t/ 


only Landau 


1.135 


0.691 


0.886 


damping 








quadrupole 


1.317 


0.791 


1.021 


damping on 








dipole and quadrupole 


1.313 


0.792 


1.020 


damping on 









Since coherent modes were excited in the experiment, 
only one bunch was simulated and compared with one 
bunch of the measured bunch signals at h = 8. The 
macro-particles of the simulated bunch were initialized 
randomly with a Gaussian probability distribution. The 
matched bunch was injected at the gap voltage of 5 kV. 
The gap voltage was modeled according to eqn. ^ where 
Atp gap is the output of the dipole feedback loop and e of 
the quadrupole feedback loop, respectively. The dipole 
feedback loop comprised an FIR filter and the dynam- 
ics of the DDS module which was modeled as an inte- 
grator (see Fig. [S] for the control loop topology). The 
input of the dipole feedback was the beam signal phase. 
The quadrupole feedback loop consisted of an FIR filter 
and an additional integral controller. The input of the 
quadrupole feedback was the absolute beam signal am- 
plitude Ai(t). The control loop parameters are given in 
Table El 

The cavity and its control loop dynamics were taken 
into account by a time constant of T cav = 20 fis. In the 
experiment, the voltage step from 5 kV to 10 kV also 
causes a small phase shift due to the nonlinearity of the 
cavity system (detuning effect). Therefore, this phase 
shift was also added to the simulation. 



E. Simulation Results and Comparison with 
Measurement 



In order to simulate the scenarios presented in section 
IIVCI using the procedure described in section IIVD1 it 
was necessary to determine the beam parameters. For all 
three cases (only Landau damping without control loops, 
only quadrupole damping, quadrupole and dipole damp- 
ing) , the rms width of the beam signal before the voltage 
jump was determined by matching a Gaussian distribu- 
tion to those measurement results which were valid before 
the voltage jump occurs. This fitting procedure leads 
to the parameters shown in Table |V] Since the bunch 
was matched before the voltage jump occurred, we have 
v x = v xQ and v y = v y0 . 

It is obvious that almost the same parameters are ob- 
tained for the three measurements which means that the 
beam quality was reproducible during the beam experi- 
ment. 

The voltage step from 5 kV to 10 kV leads to an in- 
crease of the bucket area by a factor of V2. Therefore, the 



relative bunch area vq will be 1/ \/2 times smaller than be- 
fore the step. The (x, y) phase space was defined in such 
a way that the trajectories are circles. Therefore, v x q and 
v y o will also be 1/V2 times smaller than before the step. 
In the case that no control loops are present (only Landau 
damping), this leads to v x o — 1.135/V2 = 0.803. This 
value which is valid for the new larger bucket immedi- 
ately after the voltage step was therefore approximately 
assumed in Table fl] Based on equations (121)1) and (|2"T)l . it 
leads to fs,eff = 0.8/s = 2650 Hz which is in agreement 
with the measurement. 

Now that the initial conditions for the beam were fixed, 
further parameters shown in Table HVl had to be adopted 
from the beam experiment. The only two parameters 
that were now still open were the proportional gain fac- 
tors for the two control loops. As Figures [TT] and [T2l 
show, the proper choice of these gain factors leads to a 
good agreement between measurement and simulation. 



Geometrical Interpretation for Matched 
Bunches 



For the first case (only Landau damping) we found 
v x0 = 1.135 before the voltage step occurs. According to 
equation ([27]) this leads to 

X ef f = Aip eff = 2.13. 

The (x, y) phase space leads to a bucket with the bucket 
length 2x max where x max — 7r and with the bucket height 
2y max where y m ax = 2. The bucket area equals 16. If x 
is increased from to 7r, the ratio 



y ^/2(1 — cos x) 

x x 

will decrease from 1 (circular trajectory) to 2/7T (sep- 
aratrix). If we apply this formula for x e ff, we obtain 
y e ff/xeff = 0.8214. Due to equation ([2T]) this leads to 
Vyo/v x o = 0.67. This is close to the ratio v y o/v x o = 0.61 
derived from the data in Table |Vl This is an indication 
that x e ff — 2^/v x q and y e ff — 2y/v y a may actually be 
interpreted as the effective half axes of the bunch. In our 
case this leads to a bucket filling factor of 



KXeffVeff _ TT 



16 



16 



yeff 

'eff 



= 73% 



This is close to the value 7rw /4 = 70 % based on the 
values in Table IVl 

Instead of the fitting procedure applied above, the 
value Xeff = Aifeff may also be found by determin- 
ing the two-sigma length of the bunch (in our case, the 
one-sigma length derived from the measured bunch pro- 
file is about 0.1 /is which corresponds to Aip a — 1.08). 
Based on this value, one may determine v x0 = Aip^, and 
v y o/v x0 = {y e ff I Xeff) 2 . The ratio fs, e ff/fs results from 
eqn. (1251) . Hence, all relevant parameters required for the 
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simulation can easily be determined based on the bunch 
length. 

The accuracy of the formulas presented in the paper at 
hand could still be improved by modifying the theoretical 
factors, but it is sufficient for the basic understanding and 
for the control system design. 



V. CONCLUSION 

Several specific models for longitudinal beam oscilla- 
tions have been analyzed by different methods leading to 
the following observations: 

• The definitions of the longitudinal modes in phase 
space by the bunch shape or in the frequency do- 
main by the spectral lines are not strictly equivalent 
in general. 

• In a strict sense, the term 'mode' usually implies 
linearity of the system. The examples in this paper 
show that both the dipole mode and the quadrupole 
may occur in a linear bucket whereas the sextupole 
mode requires a nonlinear bucket. 

• It was shown that the mode m = 1 shows pri- 
marily a phase modulation with the frequency fs 
whereas the mode m = 2 shows primarily an am- 
plitude modulation with the frequency 2f$. Ac- 
cording to the ODE solutions, it is also possible 
to damp these oscillations using the same type of 
modulation. Due to equation (O, a gap voltage 
modulation additionally acts in the same way as a 
phase modulation of the gap voltage whereas the 
reverse is not true. For higher orders m > 2, we 
do not find phase or amplitude modulations in the 
linear model, and it is also impossible to excite a 
sextupole oscillation. Therefore, the case m > 2 
differs significantly from the case m < 2 and must 
be analyzed using nonlinear models. 

• For control loop analysis, it is possible to use the 
quantities x and v x instead of phase and amplitude 
information. These quantities are easier to deter- 
mine since no projection of phase space onto the 
time axis is required. 

• The equations for x and v x allow an estimation of 
the area of stability for the dipole and quadrupole 
damping systems. This is not restricted to the FIR 
filters used here and could also be applied to other 
control algorithms. 

• The connection between the longitudinal coherent 
modes, the bunch mean and rms values, and the 
beam signal phase and amplitude signals has been 
demonstrated. 

• There is reason to assume that the moment ap- 
proach can be extended to the sextupole and higher 



order modes thus enabling the controller design for 
the higher order modes. 

• The commonly used transfer functions for the 
dipole and quadrupole oscillation are only approx- 
imations in the scope of our model. 

As a conclusion, the models described in the paper at 
hand allow a better understanding of longitudinal modes 
of oscillation and their damping by active feedback sys- 
tems. 



Appendix A: Linearization 

The derived equations dTUJ), CP, (HZ]), and 

x = -u s y 

tj = uj s (1 + e)(x - A<p gap ) 
v x = 2ujgV y — 2ujg(l + e)v x 
v y = -(1 + e)v x , 

can be written as a state-space model 

x = f(x,e,tp gap ) 
with the state vector 

x= (xi x 2 x 3 x 4 x 5 ) T = (x y v x v x Vy) 1 
and the nonlinear function 



f(x,e, Aip 



u>s(l + e)(xi - Aip gap ) 



gap I 



X4 

2uj s x 5 — 2w|(l + e)x 3 
-(l + e)x 4 



In the following, a linearization with Ax = x— x op around 
the operation point 



(o o vo o VQ y 



0, A(^ fl 







is performed, which corresponds to the matched circle- 
shaped bunch. This linearization (cf. Q) leads to the 
linear system 



Ax(i) = AAx(t) + bie(t) + h 2 A(p gap {t) 
with the system matrix 



(Al) 



ax 



Op 



/0 -u s 

u) S 





V o o 



\ 

' 

1 



-2w| 2w| 
-10/ 



and the input matrices 



i df 

bi = — 

oe 



= (0 -2uj 2 s v 0) 



Op 
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3 4 5 

Time/ [xs 

FIG. 14. Visualization of the function g(t) for Tr = 1 /is, 
T s = 8 /xs, At = 0.3 ms 



and 



df 



(0 -w s 0)' 



Please note that A has a block diagonal structure with 
one block corresponding to the dynamics of the bunch 
center x and one to the dynamics of the bunch variance 
v x . In addition, the bunch center is only influenced by 
<A<pgap and the bunch variance only by e. 

Comparing the equations for Ax% and Ax§ in ()A1[) 
yields Ax 3 (i) = — Ai 5 (t) and thus 

Ax 3 (t) + Ax 5 (t) = v x (t) + v y (t) - 2v = const. (A2) 

which implies that the bunch variances are connected by 
an algebraic equation and cannot be controlled indepen- 
dently. It must in principle be possible that the solution 
x of the differential equation reaches the operation point 
x op (e.g. as an initial condition). For the operation point, 



Ax 3 + Ax 5 = 



(A3) 



We define one period of a special periodic function by 

Is 1 



m = 



tT s ~ i 2 



2T, 



R 



E 

fc=i 



ak(t)a(ak(t)) - —ak(Ts) 



in the range < t < Ts {git) = elsewhere) such that 
the periodic function becomes 



-f OG 



f (0) (t)= £ 9(t-mT s ). 



a(t) denotes the unit step function. An example for g(t) 
is shown in Fig. [TJ] This function f^°'(t) is a continuous 
even function (note that g(0) — g{Ts) — 0). According 
to 



,0) 



Ts Jo 



f ip \t) cos(nujst) dt 



b^ = ^- f {p) (t) sin{nu s t)dt 
1 S Jo 

we may derive the Fourier coefficients in order to get the 
series representation 

» 



f(p)(t) = ^-+J2a^cos(nu;st)+b^sin(nu>st). (Bl) 

2 n=l 

A lengthy but straightforward calculation leads to 
2 



ai 0) = - 



T s (nia s ) 2 



Tr 



1 + £ cos (nwsTk) 



k=l 



for n^0. 



is valid which therefore holds in general for any t due to 
equation (|A2|) . 

With (|A1[) and (|A3[) . linear differential equations of 
second order can be derived for the bunch center using 
Ail = x and for the bunch variance using Ax% = v x — vo- 



x + uj s x = ojgAtp 



.</«/< 



v x + Aw 2 s {v x - vo) = -2u 2 s voe. 



Appendix B: Fourier Series of the Dipole Oscillation 

Signal 

As an abbreviation, we set 



T,, = AT,,. + A I .si,, { 2;rfcp- ) Al ■: T„i2 



and 



ak(t) = t-T k - 



,(0) 



r 2 
6T R 



Is. 



E \ T,: 

k=l ~ 



' k - T k ) forn = 0, 



= for all 
It is obvious that for n ^ 0, 

k 0) i< 



TR.(nu} S y 



holds. Hence, we have |a£°^| < Mn K with real con- 
stants M, k. According to [13], Corollary 2.4-3b, the se- 
ries f( p '(t) in equation (|Blj) converges in the space V of 
distributions. In V it is allowed to differentiate this se- 
ries without further requirements ([lOj, Corollary 2.4-3a). 
For the function g(t), we obtain 



g'(t) 



T s ~2t 
2T R 



E 

k=l 



a{fl k {t)) - —a k (T s ) 
J-S 
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(for < t < T s , g'(t) = elsewhere). We find 



with 



T R 



g'(T s -) 



y - 



k=i 



such that 



'2T R 



E 

fe=i 



7\ 



f (1) (t)= 9'(t-mT s ) 



m— — oo 



(regarded as a locally integrable function) has unit steps 
at t = mTg. Both g' and are locally integrable 
functions and may therefore be regarded as regular dis- 
tributions. 

As mentioned above, we may now differentiate this dis- 
tribution. Thus we obtain in V: 

f (2) (t) = -—+ y 6 M)) 

K— — 00 

On the other hand, the derivative of cqn. fjBip yields 

oo 

n=l 
oo 

= — y^ a^\nuj s) 2 cos{nuj si) 



,(2) _ ±_ 



i + y^ c ° s ( nu, sTk 



k=l 



for n^O. 



This completes the proof of eqn. 
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